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Abstract 

We calculate electromagnetic production from highly excited hadronic 
matter created in heavy ion collisions. The rates include the usual lowest or- 
der processes in quark-gluon plasma plus the usual reactions in the hadronic 
phase, related with p, a± mesons. The space-time integration is done us- 
ing a hydrodynamical model. We have found that dilepton mass spectrum 
agrees with results of other previous works, but disagrees with the CERES 
dilepton data. In order to explain these data, some "unconventional" pro- 
duction mechanism need to be incorporated: we discuss especially the notion 
of modified p, ai masses, which indeed may explain the data. Other sugges- 
tions (e.g. the longer-lived fireball) to increase the production of low mass 
dileptons seem to be insufficient for the task. The results for direct photon 
production are below the current WA80 experimental bounds, for all variants 
considered. 
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1 Introduction 



The main goal of heavy- ion collision program in the AGS/SPS energy range 
(10-200 GeV/A) is to produce hot/dense hadronic matter with the energy 
density of the order of few GeV/ fm 3 and to study its properties. Especially 
interesting are early stages of the collisions, when theory predicts existence of 
the QCD phase transition into a new phase, called the Quark-Gluon Plasma 
(QGP). However, so far no direct experimental evidence of the QGP has been 
found. The main reason for this is well known: strong collective interaction 
in the system as it expands and cools erases most of the traces of the dense 
stage. As a result, the observed hadrons come mostly from a dilute freeze-out 
stage, with a final temperature Tf = 120 — 140MeV.[]I] 

One possible way to study the earlier stages (to be discussed in this paper) 
is to look for phenomena which mostly happen very soon after the collision, 
such as production of dileptons and p/ioton,s00 . At high (RHIC/LHC) en- 
ergies one may hope for some kinematic enhancement of the QGP signals 
because the initial stage is much hotter than hadronic mater, T{ » T c 
H £|, but this is certainly not the case for AGS/SPS energies. Therefore, in 
this energy domain the main signals for the phase transition still come from 
the hadronic stage, not the QGP. 

Experiments designed to observe the direct photons or dilepton continuum 
produced by excited hadronic matter in heavy ion collisions are generally 
much more difficult to perform compared with measurements of hadronic 
observables. Therefore, only recently were the first photon and dilepton 
measurements announced by four CERN SPS experiments: NA34/3, NA38 

3 Another possibility is to look for signals which are accumulated during the evolu- 
tion: the well known examples include excessive production of strangeness or charmonium 
suppression. 
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and CERES (NA45) for dileptons and WA80 for photons. It was found 
that dilepton production exceeds backgrounds expected from hadronic and 
charm decays. Furthermore, the signal also exceeds theoretical expectations 
for "conventional" processes, both in hadronic and quark-gluon matter ||. 
Especially dramatic is the excess observed by CERES || in the mass region 
M e + e - = 0.3 — 0.6 GeV. This observation has since created a rapidly growing 
theoretical literature. 



Calculation of the dilepton/photon yield consists of two components: (i) 
evaluation of production rates (see section |2]); and (ii) their integration over 
the space-time evolution of the collision (see section |3|). In the so called 
conservative approach (the well-known hadronic processes with vacuum pa- 
rameters and the usual space-time evolution of heavy-ion collisions) several 
groups have obtained rather similar results, which however do not explain 
the CERES data, neither in magnitude nor even in the shape of the mass 
spectrum. 



This situation has lead to many "unconventional" hypotheses, which in- 
clude (in a more or less chronological order): (i) dropping m p ||, |7|; (ii) high 
pion occupation numbers at low momenta ||; (iii) a very long-lived fireball 
[|5[]; (iv) dropping m v > [M, [TTL |T2[; (v) a modified pion dispersion curve JTB 
(vi) dropping m ai (discussed below). 



Let us start with the possibility (ii), which is usually taken into account by 
introduction of the pion chemical potential which is approximately equal 
to m n . Although the true nature of the low-p t pions is not yet completely 
clear, most probably they come from the resonance decay and/or spectral 
modification due to collective potentials. Both are late-stage phenomena, 
which can hardly affect the early-stage dilepton production. Furthermore, 
studies of this explanation made in JF| have shown, that the low-M dilepton 
enhancement due to //^ ~ is way too small compared to CERES data. 



In section || we will look at (i), the T-dependent m p , and (similar to 
conclude that it may indeed describe the data. Furthermore, we have linked 
it to (vi) by theoretical arguments related to chiral symmetry restoration and 
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have derived experimental consequences of the most plausible scenario, both 
for dilepton and photon (section || ) yields. 



we 



We also study separately a proposal (iii). In our previous paper JT3 
found that (at least in a hydro approach) the so called "softest point" of the 
Equation of State (EOS) leads to especially long-lived fireball. Although it 
is expected to happen at collision energies way below those for CERN exper- 
iments, in section [7] we have pushed this idea to the extreme and assumed 
the scenario with the long-lived fireball. The results for the dileptons are 
quite disappointing: very different scenarios of space-time evolutions give 
very similar dilepton mass spectra. So, with the standard dilepton rates and 
fixed masses it is not possible to explain the observed dilepton yield by longer 
lifetime. 



2 Electromagnetic processes in hadronic mat- 
ter 



Dilepton/photon production in the QGP phase is based on fundamental QCD 
processes like qq — > e + e~ H and was calculated long ago. The rate in the 



pion gas due to tttt annihilation was considered in [15], and those two basic 



processes can be included by the "standard rate" formula: 



dR a 2 



d 4 q 487T' 



Fe~T (1) 



where the rate R is counted per unit volume per unit time, q is 4-momentum 
of the virtual photon (q 2 = M 2 +e ~ = M 2 ), F is a constant in QGP and the 
usual pion form-factor in the pion gas, which can be written in standard 



4 



vector- dominance form [] 



F h = (Hadronic) 

I (2) 

def n ^ of-, , 2ml \ ( 4m^ \ 2 



12V eMl + ^ (QGP) 



Later additional processes including p mesons were added. It was also pointed 
out in |16| that a\ meson is very important^, especially for photons and low- 
mass dileptons. Further work was done in refs. |T7|, |I5f . 

Before we come to specific formulae, let us remark on some misunder- 
standing of the role of the resonances, which has even led to double counting 
in some previous papers. Unstable particles like p meson can be considered 
either (i) as parents of decay processes ( such as p — > e + e~ , a\ — > e + e~n ) or 
(ii) as intermediate states in particular reactions with "stable" particles (e.g. 
7r + 7r~ — > e + e~). Furthermore, the same resonance can enter as an interme- 
diate stage in many different reactions: This created the impression that by 
considering many of these reactions, one can in fact increase the dilepton 
yield. 

However, this is not true. In thermal equilibrium the average number of 
mesons depends on their mass but should not depend on their width (which 
only shows how often these particles are created and decay). In order to 
understand how this happens (we'll give detailed discussion in the Appendix), 
recall that standard Breit-Wigner amplitude is proportional to T in T out /((M— 
M res ) 2 + T^/4). Summing over all possible "in" channel one gets T tot in 
the numerator, and after that one may approximately^ substitute the Breit- 
Wigner amplitude simply by T out 5(M — M res ). The resulting rate is nothing 

4 For the quark masses, we adopt the perturbative result |l5| m q (T) — -^=, with g = 2.05 

which corresponds to a s — 0.33. 

5 In order to explain why a\ is important, let us go "backward in time": it is the first 
hadronic resonance which may be excited in a collision of a photon, real or virtual, with 
a pion. 

6 Provided the total width is not large compared with the temperature! 
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else but the decay contribution (i) mentioned above: clearly one should not 
include it twice. We have found it rather convenient for our applications to 



Dilepton Rate 



l shil'W. 1 = 160 MeV 



Dilepton Rate 




Figure 1: Dilepton production rates (a) standard (b) comparison with max- 
imally shifted according to variant D discussed in section [j| 



use expressions recently derived in [19|: They contain processes in zeroth and 



first order in the pion density. The expression for the dilepton rate reads 



dR 
~¥q 



a 



3n 3 q 2 1 + exp(q /T) 



3g 2 imIT(g 5 



d 3 kn(k ) 



(2tt)32A; /2 



[-12q 2 ImU v (q 2 ) + Q(k + q) 2 ImU a ((k + q) : 
+ Q(k-q) 2 ImTl a ((k-q) 2 )] 



(3) 



where k is the pion and q is the virtual photon momentum. Here ImU v , ImU a 
are imaginary parts (or spectral densities) for vector and axial currents. If 
this expression is understood as expansion in the pion density, they should 
be evaluated in vacuum, and thus related to experimental data on e + e _ — ► 
hadrons and r lepton decay. Furthermore, ImU v , ImU a can be approxi- 
mated by well-known contributions of p, a\ resonances, which produce dilep- 
tons by their decays into e + e~ and 7re + e~, respectively. We chose the follow- 
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ing parametrizations for ImU v and ImU a (see section |5|) 



Imli v 



m 2 p0 m p T p 



ImU-a 



M 2 (M 2 - m 2 ) 2 + m 2 p T 2 p 
fgm a T a 

■ ™l) 2 + ™ 2 cF 2 a 



(M 2 



(4) 
(5) 



where f w =93 MeV, T f 



Mev, r a = r 



«0" 



M 2 



aO 



r p0 ^, r> = 149 MeV, m p0 = 770 MeV, f a =190 

pO 

: 400 MeV, m a0 = 1210 MeV g| • 



The resulting rates |T9[ are shown in Fig.|T](a), where we show the con- 
tributions of the vector and axial parts separately. We have used simple 
Breit-Wigner parametrizationT|with vacuum parameters for resonances taken 
from Particle Data Tables, while part (b) corresponds to both p, a\ masses 
shifted to \m p o in the "mixed phase" . (We assume the critical temperature 
T c =160 MeV, and the mass shifts will be discussed in detail below). In both 
cases, the direct channel resonance p dominates around its mass, while the 
ai contribution takes over at small masses. 



3 The model for the space-time evolution 



The hydrodynamical model was suggested by Landau more than 40 years ago, 
and its applications to heavy ion collisions has a long history. Experiments 
with heavy ions at lower energies (~ 1 GeV/A) were able to detect collec- 
tive motion of nuclear matter by comparing velocity distribution of different 
nuclear fragments. 

One important finding is that in the AGS/SPS energy range for heaviest 
nuclei (Au Au at AGS and Pb Pb at SPS) the rapidity spectra of n,K,p,d 

7 This is why vector contributions does not vanish below 2m v . 
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are consistent [[IJ with a simple hydro description: a convolution of thermal 
motion at breakup (which depends on the particle mass) with a collective 
flow common for all species including baryons^. An open question however is 
whether thermalization is rapid enough, so that one can use hydro description 
from very early times, or collective motion is formed only at later stages. In 
this work we assume that the former is the case, and will use hydro description 
for dilepton/photon production. 



The second important observation is related with the equation of state 
(EOS). For finite T and zero baryon density, a conventional parametriza- 
tion of the EOS can be given by a "resonance gas" [20| below T c , plus a 
bag-model QGP above T c . We fixed the phase transition point (smoothened 
for numerical purposes) at T c = 160 MeV. In hadronic phase the speed of 
sound dp/de = c 2 s = 0.19, and in the plasma phase the bag constant is 
B = 0.32 GeV/fm 3 . The results are plotted in Fig.§. This EOS is in rea- 
sonable agreement with the lattice results. Unfortunately, there is so far 
little progress in lattice simulations for non-zero baryon density. People have 
extrapolated various models which are successful at nuclear densities, like 
Walecka model, but it is unclear how well such extrapolation works. Some 
guidance can probably be provided by hadronic gas including baryonic reso- 
nances, for example the one discussed in U . In Fig.0 we have shown the cor- 
responding curve (dashed line) for baryon chemical potential = .54 GeV, 
corresponding to AGS breakup conditions. In this case the baryon/meson 
ratio is about 1, one order of magnitude larger than baryon admixture in 
central region at SPS energies: and still the EOS in p, e coordinates looks 
approximately the same as the one we adopted. 



In summary, the first observation shows that for heavy ions the flow of 
the entropy (mesons) and the baryonic charge (nucleons) are about the same. 
The second shows that EOS (in p(e) form) is not much affected by baryonic 
charge. Taken together, they provide a reason to ignore baryonic charge in 
our hydro calculations. 



8 In the framework of cascade models such as RQMD this topic was studied (2TJ and it 
was also concluded that if one cut the excited system into elements, the mean velocity of 
different species are the same with accuracy 10-20 percent. 
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Equation of State 



0.20 




Figure 2: The EOS in hydro-relevant coordinates: the ratio of pressure to 
the energy density p/e versus e. Notice the minimum, which we refer to as 
the "softest point". The dashed line is the EOS for finite baryon chemical 
potential //& = .54 GeV from Braun-Munzinger and Stachel. 

We use standard equations of (non- viscous) relativistic hydrodynamics for 
central collisions (i.e. we assume axial symmetry). The hydro equations were 
solved numerically using the first-order Lax finite difference scheme. Energy 
and entropy conservation is monitored, and we have also made comparison 
with results of several earlier works was made to ensure that technical aspects 
are under control. 

The major uncertainties one faces dealing with hydrodynamical mod- 
els are the initial conditions^. Those can be described by the following set 
of parameters: (i) the initial size z ; (ii) the fraction of thermalized en- 
ergy/collision energy k; and (iii) v ijZ describing initial distribution of the 



9 Although cascade-type event generators (Venus, RQMD or ARC) provide some guid- 
ance, their physical basis is questionable exactly at the first 1-2 frn/c of the collision, when 
we need it. 
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longitudinal velocity 



v z (t = 0) = Vi >z * tanh(z/zo) (6) 



We consider specifically two cases of central collisions: (i) S-Au 200 GeV/N 
and (ii) Pb Au 160 GeV/N for which the dilepton/photon data were taken. 
In all cases considered, at t = we assume that some part k of total energy 
goes into the thermalized matter at rest, and fix k demanding that at the 
end of the expansion (at Tf = 140MeV) the predicted number of pions0 is 
the same as observed. 

One extreme scenario is complete equilibration, Vi >z = a la Fermi- 
Landau. Furthermore, if thermalization is very rapid, the usual Lorentz 
contraction is enhanced by the ordinary compression of matter by shocks, 
leading to very high e» ~ 10 GeV/fm 3 initial energy density, well above the 



phase transition region. As shown e.g. in [£3[ this scenario is incompatible 
with SPS data on rapidity distribution. However, as shown in Fig.|3](a), for 
Vi.z = 0.9 one obtains reasonable description of the longitudinal motion^. 



4 Conventional dilepton yields 



Now we evaluate the dilepton yield, by integrating the production rate over 
the space-time. The integral is taken between the surface at which initial con- 
ditions are set (the so called "pre-equilibrium" contribution is thus left over) 
and the "breakup" surface, at which the density is so small that secondaries 



10 Multiple studies have shown that about 1/3 of pions come from resonance decays. 
We included this fact in the normalization. 

11 Distribution of other secondaries, as well as the transverse motion we plan to present 



elsewhere. The hardon rapidity data is taken from 1 24 (The data actually refers to negative 

dN 
dy 



hadrons. We scaled them by 0.7 to account for resonance decays. The y > 2.65 data 



for S+Au has been reflected about y = 2.65 to give data for y < 2.65.). 
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Figure 3: Typical hydro evolutions for 200A GeV S+Au (upper part) and 
160A GeV Pb+Au (lower part). The left part shows hydrodynamical solu- 
tion in the plane time- longitudinal coordinate. Solid lines are lines of constant 
energy density, the dotted ones correspond to constant longitudinal velocity. 
Thicker lines show the end of the mixed phase (e = 1.47GeV/fm 3 ) and the 
break-up conditions (e = 0.31GeV/fm 3 ): at z = that happens at about 
time 5 and 7 fm/c, respectively. The figures on the right shows the calcu- 
lated pion rapidity spectra (lines) compared with data, where the data has 
been scaled by 0.7 to account approximately for resonances. 

can fly away without interaction 0. The surface should be determined from 

12 In principle, some pion annihilation can happen even later: we have evaluated this 
contribution and have found that it is small. Note also that p, uj, rj etc decays after breakup 
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the mean free path, and for pions it corresponds to Tb rea k up = 140MeV r . Fur- 
thermore, in order to compare with CERES data one has to apply their exper- 
imental cuts: For S+Au, p T > 200 MeV/c; 6 ee > 35 mrad; 2.1 < r] < 2.65. 
For Pb+Au, p T > 175 MeV/c; Q ee > 35 mrad; 2.1 < r, < 2.65. 

The results we obtained are shown in Fig]4|(a): they are split into 4 con- 
tributions: from pure QGP phase (denoted by Q); QGP part of the mixed 
phase (Mix/Q); hadronic part of the mixed phase (Mix/H); and finally from 
the hadronic phase (H). Note that, as expected, QGP contribution domi- 
nates at high masses, while hadronic contribution dominates at p region and 
below. Although in the rates shown in Fig. |T] one can see quite substantial 
contributions at low masses from a\ decay, it is reduced by the small exper- 
imental acceptance in this region so dramatically that it falls well below the 
data. 

For check of consistency, we have compared our results with those of other 
approaches in Fig.f|(b). All three use "conventional rates" which are similar 
to ours (except the a\ part) but very different dynamics of the evolution. Li- 
Ko-Brown approach is based on cascades started from RQMD-based initial 
conditions^. Srivastava et.al. assumed Bj or ken-scaling longitudinal hydro- 
dynamics, with solved hydro in the transverse planeQ. We conclude that the 
agreement between 3 papers is reasonable in terms of the total yield, and 
even very good in terms of the shape of the mass spectrum. Since all three 
results disagree with data, one may conclude that it is not possible to explain 
CERES data by "conventional sources" . 

In Fig.[|(c) we also show our predictions for 160A GeV Pb+Au collisions. 
We compare them with new preliminary CERES data reported at Quark 
Matter 96 ||. Since the data are for non-central collisions, we have multi- 
plied our curve by the factor 1.7, the ratio of the observed dilepton yield for 
minimum bias sample < n c h >— 260 to that in central sample < n c h >~ 400. 

are precisely the "hadronic background" , which was separately calculated, by the experi- 
mental group itself. 

13 Note that it includes additionally the <j) contribution ignored by us and others. 
14 Note that this curve is already corrected for the double counting which was present 
in the original paper, but their results are still somewhat larger. 
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In this case the disagreement is not statistically as significant as in the pre- 
vious reaction, but still it seems that the measured shape of the spectrum is 
different from the predicted one. 




.5 1 1.5 o 05 1 

M (GeV/c ) M <°= v > 



Figure 4: CERES dilepton spectra from conventional sources. (Top) our 
conventional spectra for 200A GeV S+Au, shown without and with back- 
ground from hadronic decays; (Bottom) our results for 160 A GeV Pb+Au 
w/o and with background. 
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5 Dilepton yield from modified /?, a\ mesons 



> 



o 
o 



-a 



Case A 
Case B 
Case C 
Case D 
Case E 




M (GeV) 

Figure 5: Dilepton yield in 200A GeV S+Au collision for different possible 
scenario of chiral restoration in m p — m ai plane. 

Dileptons, unlike secondary hadrons, are produced at relatively early 
stages of the collisions One possible explanation of why the "standard" rates 
fail to reproduce the observed excess of low-mass dileptons may be a long- 
debated idea: hadronic properties may be modified in high density hadronic 
matter. 
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Figure 6: Direct dilepton yields (without and with background) for 200 A 
GeV S+Au (top) and 160A GeV Pb+Au (bottom) using mass scheme D. 



We know definitely that the nucleon mass is modified in nuclear matter, 
and shifts of the vector meson masses in it are debated for a long time, see 
e.g. 



251 EB| 



In a low-T hadronic gas (which may include an admixture of baryons) 
one can relate modification of mesons (e.g. of p) to the up (and Np) for- 
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ward scattering amplitudes [E7J . This approach predicts certain momentum- 
dependent optical potential, which can be loosely interpreted as (relatively 
modest) shift of m p downward. Finite T/density QCD sum rules (see e.g.pSfl 
and references therein) relate hadronic properties to the quark condensate 
< qq >, which decreases toward chiral restoration transition. This reasoning 
has culminated in the so called Brown-Rho scaling idea, according to which 
all hadronic dimensional quantities get their scale from < qq >: therefore all 
masses are predicted to vanish at T — > T c . 



In the instanton model the chiral restoration is due to transition from 
random instanton liquid to a gas of instanton- anti- instanton molecules [EU 



the latter survive the phase transition and lead to new type of quark inter- 
actions unrelated to < qq >. Although results of available simulations PC 



have not reported any definite conclusions about the rho meson mass, at 
quark level it was clearly demonstrated that at T > T c an effective quark 
mass is substituted by "effective energy" (or "chiral mass") of comparable 
magnitude, which imply that hadronic masses should not vanish at T —> T c . 



The matter however is by no means settled, and quite different sugges- 
tions about hadronic properties close to (or even above) T c can be found in 
literature. For example, effective Lagrangians lead to a prediction of rising 
m p (T) PHI , moving it about half way toward the mass of its chiral partner 
a±. Furthermore, at SPS about 1/10 of secondaries are baryons, and thus 
finite baryon density effects should be added to the finite-T shifts mentioned 
above. 



Therefore, we have decided to proceed empirically, assuming various sce- 
narios of hadronic mass evolution in dense matter. A particular point we 
want to make in this paper is that "dropping m p (T)" idea is consistent with 
chiral symmetry restoration only if appropriate modifications of its axial 
partner m ai follow. Strict theoretical relation between the two were derived 



via Weinberg- type sum rules [32]. Possible scenarios of how chiral restora- 
tion may proceed are therefore shown in the m p — m ai plane in Fig||(a). 
For example, path B corresponds to results of Pisarski |31| , while E corre- 



sponds to Brown-Rho scaling. Anyway, both mesons should become iden- 
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tical at T c , so the path of thermal evolution should end up at the diagonal 

m p (T c ) = m ai (T c ). 

To test these ideas, we need a simple consistent model of how a changing 
p mass can change the dilepton production rates discussed in section ^|. In 
this respect, the standard expression for the pion form-factor (§) happen to 
be quite misleading: it includes m p both in numerator and in denominator, 
but simply to make it T-dependent everywhere would in fact be a mistake. 
Recall the usual reasoning why there is m A p in numerator is that at M = 
the form-factor should be F(0) = 1. This statement should hold for the the p 
contribution only if one assumes the so called vector dominance, demanding 
that the whole form-factor (and not just a part of it) is given by the rho 
pole. It is quite accurately satisfied at T=0, but for T-dependent m p vector 
dominance has no reason to persist! 

Actually the numerator of the form- factor contains a combination m 2 p Y ' e + e -T p , 
and in what follows we assume that both widths are 7 "-independent. Note 
that only assumption T e + e - = const{T) is actually important, because (i) 
m 2 cancels when one returns to non-relativistic form of Breit-Wigner; and 
(ii) as we have argued above, the value of the V p is nearly irrelevant for the 
rate anyway. 

Let us concentrate only on the p part of the rate and write the dilepton 
production rate as (see appendix): 



dR 




+ (M)' 



(7) 



where m 2 = m 2 pQ + i?ell and II is the p self-energy. One can define the 
M-dependent width of p by the relation ImliiT = 0) = f — m P0 T(M) which 
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is generalized to non-zero T by ImU(T, M) = —m p (T, M)T(M) where we 
ignored the modification of the width with T. Finally we can rewrite Eq. (^) 

as 



dR ^fiKfi+M. 

d 4 k ~ 7T 3 M 2 \ M 2 

m p {T)Y{M) 



(M 2 - m 2 f + m 2 p (T)T 2 (M) Ve^ - 1 




(8) 



where we have suppressed the M-dependence of m p , m p and we used the 
relation f 2 = ^p-. 

For simplicity we identify m p (t) and m p (T), with both given by 
m{T) a 



m T/T c 



+ r '^~, ( 9 ) 



where r is the value of when T = T c , a characterizes the abruptness of 
the mass shift, while b is chosen to make = 1 at T = Tt = 140 MeV. 

' mo J 

Now we can test whether a T-dependent masses can indeed describe the 
CERES data in the lower mass region. Integrating over space-time according 
to hydro calculations described above we get the dilepton mass spectra as 
shown in Fig. |5|. The variant D, with m p (T c ) = m ai (T c ) ~ |m p (0), does the 
best job, and Fig.^| we show contribution of separate stages in this scenario, 
for both 200 A GeV S+Au and 160A GeV Pb+Au. It is the hadronic part of 
the mixed phase which is responsible for the observed excess at M = 0.2 — 0.6 
GeV. Furthermore, depending on how exactly m p (T) goes to its limit at T c , 
one can change the shape of the resulting mass spectrum. This statement is 
demonstrated in Fig.[7|, where (a) shows several scenarios, with the resulting 
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mass spectra in Fig.|7|(b). The case E in ^ which corresponds to massless 
p and di, generates more low-mass pairs than the hadronic cocktail. Better 
data of r] production is needed to rule out this case. 

Of course, specific dynamical models lead to more complicated picture: 
instead of simple universal m p {T) for any p meson, they lead to specific 
modification of the whole dispersion curve. Clearly, mesons which travel 
fast relative to matter are modified differently from those which have zero 
velocity. Furthermore, since the meson-meson and meson-baryon scattering 
is dominated by resonances, this dependence may even be non-monotonous. 
Future high-statistics studies may look into those matters by considering 
mass spectrum at different rapidity and/or p t of the dileptons. 

In conclusion: we have demonstrated that the CERES data can be fitted 
pretty well over the whole experimental mass range by assuming a particu- 
lar m p (T), without any change in the production mechanism or space-time 
evolution. 




Figure 7: Different temperature profiles for m p (T) (left) The corresponding 
dilepton spectra, (right) For meaning of the parameters a and r, refer to 
equation (9) 
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6 Searching for the modified a\ contribution 



Great interest related with modified hadronic masses is explained by a possi- 
ble relation between this phenomenon and the chiral phase transition. As it 
was repeatedly emphasized above, modification of p imply certain modifica- 
tions of cii as well. In this section we discuss how this can be experimentally 
verified as well. 



The important role of a\ for production of photons and dileptons was 
discussed in |16|, |], |17[ [18[] As we have shown in section |], the main a± 
contribution is at low dilepton masses. However, CERES has adopted cuts 
which makes its acceptance very small at small M e + e -, and therefore, the a\ 
contribution to its mass spectrum is not very significant. Those cuts were 
made in order to get rid of hadronic background, such as rj Dalitz decay. So, 
by simply looking at small mass region one cannot find the a x contribution. 



Nevertheless, one can try to locate a kinematical "window" in which it 
may be better seen. We have found that the best range of invariant dilepton 
masses is M e + e - fa 300MeV. Above it hadronic background is relatively 
small, but the contribution of p channel and even of the quark annihilation 
in the QGP phase exceed the a\ contribution. Below this mass they are not 
so important, but soon the hadronic background from ordinary Dalitz decays 
take over. 



Furthermore, there is a particular p t dependence of the effect. Consider 
first the unmodified oi, decaying into ne + e~ . For relatively small dilepton 
masses, in the a\ rest frame both the pion and dilepton get half of its total 
mass. Including thermal motion and (rather large) a\ width, one still finds 
a broad maximum in dilepton production rate at dilepton energy about 600 
MeV. Integrated over longitudinal momentum, one gets dilepton p t distribu- 
tion shown in Fig||(a), with a maximum. It looks very different compared to 
the dileptons coming from p decay and having the usual thermal p t depen- 
dence. Now comes the main point: if the ai mass shifts down, as expected 
from chiral restoration, then the wide peak is absent and the shape is differ- 
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ent: see Fig.|](b). 



In summary, it is important to locate and to test properties of a\ con- 
tribution. We have shown how sensitive dilepton production is to its modi- 
fication in dense matter. Although this task is not easy, next generation of 
experiments (with accurate knowledge of hadronic backgrounds, measured 
dependence on the event centralities etc) can in principle do it. 




Figure 8: Dilepton p t spectrum for 200A GeV S+Au. The top two figures 
are for mass scheme A (i.e. unshifted masses) while the bottom two are for 
mass scheme D. 



21 



7 Can one explain the dilepton excess by a 
long-lived fireball? 



In our previous work ]I4[ it was shown, that under certain conditions the 
"softness" of EOS near the phase transition can affect the longitudinal ex- 
pansion, so that even the global lifetime of the excited system increases sub- 
stantially. It happens in a window of collision energies such that the initial 
energy density is close to the "softest point", which makes secondary accel- 
eration of matter is impossible. Due to uncertain initial conditions, we do 
not exactly know to which collision energies this window correspond, but 
estimates put it in the region E L ab ~ 30GeV/A, between AGS and SPS. 
Some examples of the solution we got at higher energies and at the softest 
point are shown in Fig|| One can see that the long-lived fireball is slowly 
burning, with the lifetime reaching more than 30 fm/c! 



It was further been proposed that one can test these unusual predictions 
experimentally by: (i) Looking for the maximal lifetime (or the minimum of 
the HBT parameter A) ||14j| ; (ii) Looking for the minimum of the "directed" 
flow in the collision plane |33[ ; (iii) Looking for the nearly isotropic distribu- 
tion of dileptons, produced in the long-lived fireball . In connection with 
(i) it is very intriguing that the E802 AGS experiment reported preliminary 
studies of HBT which indicate significant (~ 40% ) growth of lifetime for the 
most central Au Au collisions l34| . 



We now would like to check whether the long-lived fireball, if present at 
SPS energies, can enhance the dilepton production to a degree necessary to 
explain the CERES data. For that we take the initial conditions as shown 
in Table [T], which indeed give the initial energy density close to the softest 
point value. 



The results however indicate that this scenario can neither be made con- 
sistent with rapidity distribution of secondaries, nor does it actually result in 
larger production of dileptons. In Fig.|T0| we compare two expansion scenarios 
with unmodified masses. The reason for this is as follows: longer lifetime is 
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ENERGY DENSITY £ at r = ENERGY DENSITY £ at z = 
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Figure 9: Space-time evolutions for Standard Bjorken-like scenario (top) 
Long-lived fireball scenario (bottom) 



compensated by smaller 3-volume. And also, the shape of the dilepton mass 
spectrum is different from the observed one. 

Furthermore, if masses are modified (say according to variant D of Fig.[5](a)), 



both space-time pictures become compatible with data, see Fig. 11 
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Table 1: Parameters for the two hydro models used 

Bjorken-like Long-lived Fireball 



Initial energy density Cj(GeV/fm 3 ) 


9.0 


1.5 


Initial temperature Tj(GeV) 


0.27 


0.16 


Total energy in fireball E^GeV) 


830 


930 


Initial long, velocity v ZQ /c 


0.0 


0.0 


Initial transverse vel. v ro /c 


0.0 


0.0 


Initial long, half-size zo(ftri) 


1.1 


8.0 


Initial transverse radius ro(fm) 


3.8 


3.8 


Critical temperature T c (GeV) 


0.16 


0.16 


Freeze-out temperature 7/ (GeV) 


0.14 


0.14 




Figure 10: 200 A GeV S+Au direct dilepton yields using unshifted masses for 
the usual Bjorken-like expansion (left), and for the long-lived fireball (right) 

8 Photon production 



Let us briefly comment on the history of theoretical and experimental studies 
for photon production. As for dileptons, the radiation from QGP phase was 
already calculated in 0. Reactions in a gas of tt, p mesons were considered in 
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Figure 11: 200A GeV S+Au direct dilepton mass spectrum using the p, a\ 
masses in scenario D for Bjorken-like hydro (left) and for the long-lived fire- 
ball, (right) 



fL5| , and importance of a\ — > 77? reactions was pointed out in . Following 
our dilepton calculations, we use the photon rate given in |T9 . 



Experimental observation of direct photons is an extremely difficult task. 
Unlike dileptons, we cannot tune the invariant mass and therefore for all 
photon momenta the background from hadronic decays (mostly 7r°, 77) domi- 
nates the signal of the "direct" photons. Therefore, the issue is very accurate 
measurements of these backgrounds and inclusive photon spectra, with sub- 
sequent subtraction. At the moment, only upper bounds on the direct photon 
cross section has been given by the WA80 experiment |3§ 



15 Earlier preliminary WA80 data indicated a non-zero effect: its comparison with the- 
oretical expectations based on "standard sources " p] lead to conclusion that those data 
significantly exceeded the expectations. Later in J22| this conclusion was disputed: larger 
photon yield (now consistent with data) was obtained. The issue under debate was ba- 
sically how one should normalize the hydrodynamical initial conditions. Eventually, the 
preliminary data were withdrawn after the reanalysis. 
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In Fig.[12] we show results of our calculations of the direct photon pro- 
duction. They correspond to unmodified hadronic parameters, and are per- 
formed both for standard space-time scenario (a) and the long-lived fireball 
(b). The main conclusions are: (i) The theoretical predictions are in both 
cases below the experimental bound in the whole region, although the differ- 
ence between them is not that large, (ii) The two scenarios show significant 
difference at large p t ~ 3GeV, mostly due to existence of relatively hot QGP 
in the first case. It would be very important to pursue the issue further, in 
WA80 or elsewhere, and try to observe radiation from QGP and hot hadronic 
gas. 




Figure 12: Direct photon production in 200A GeV S+Au compared with 
preliminary WA80 upper bounds on direct photon production. Our predic- 
tions are for the Bjorken-like expansion (left) with hot initial stage, and for 
the long-lived fireball, (right) 



9 Conclusions 



In summary, we have studied dilepton/photon production from highly excited 
hadronic matter produced in heavy ion collisions. 
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Using the rates from the usual hadronic reactions with vacuum param- 
eters and the usual hydro description of space-time picture, we obtained 
results which are consistent with those of previous works but cannot account 
for low-mass dilepton excess observed by CERES experiment. However, if 
the p mass is shifted to about 1/2 of its value around the critical temper- 
ature, the data can be explained. More detailed observed distribution (e.g. 
over dilepton mass and p t ) should provide better understanding of whether 
this explanation of the observed dilepton excess is in fact correct. 

Furthermore, we have pointed out that chiral restoration demands that 
the mass of a\ meson should be shifted together with the mass of p In order to 
get more convincing evidence that the excited matter is indeed approaching 
chiral restoration, experimental observation of the ai-related component is 
essential. We discussed in which kinematical window one should look for it, 
and evaluated the magnitude of the effect. 

Among suggestions to explain the dilepton excess was a proposal that the 
space-time picture of heavy ion collisions can in fact be different, with pos- 
sibly a longer-lived fireball. However, although for specific initial conditions 
this long-lived scenario is hydrodynamically possible, it does not significantly 
enhanced production of dileptons because longer lifetime is compensated by 
smaller spatial volume. Furthermore, this scenario is incompatible with the 
observed spectra of secondaries. 

Finally, all scenario considered lead to photon production well below cur- 
rent upper limit on direct photons from WA80. The results however are 
sensitive to expansion scenario. The QGP component in particular becomes 
dominant around p t ~ 3GeV. 
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Appendix: Thermal Dilepton Rate 



It can be shown |36| that the thermal rate of dilepton production in a hot 
hadron gas is: 



dR 



2e 2 



r d?p + fd 3 p. 4 
J ~E~ J ~e7 



e f3w _ 1 



(2vr) 6 (fc 2 ) 



V+V- + V+V- ~ 9 >MU (P+ ■ P- 



m, 



(10) 



Using vector dominance to relate the imaginary part of the photon self- 
energy n^j, and the imaginary part of the p propagator A^, ]36| , then inte- 
grating over p + and p- , we get 



dR 



a 2 m 4 p 



d 4 k 3tt 3 / 2 M 4 



1 + 



2m 2 
M 2 " 



(11) 



This can be manipulated further into the form [37] [36 



dR 2ujdR 



d 4 k dM 2 d 3 k 7T 3 f 2 M 2 



9 4 

-a 2 rrr p 



1 + 



x- 



ImU 



(M 2 -m 2 ) + (ImUf 



2m 2 \ 
W 2 ) 
1 

_ I 



Am 2 \ * 



M 2 



(12) 



where rh 2 = m 2 + ReU and II is now the transverse or longitudinal part of 
the p self-energy. We'd now demonstrate that the above rate can be written 
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in terms of thermal p decay rates, we are interested in the invariant mass 
distribution (assuming IT = H(M)), 



dR —a 2 mi 



dM 2 2tt 3 f 2 M 2 



1 + 



2m 2 ' 
~M 2 



i 

Am 2 \ 2 



M 2 



ImU. 



(M 2 -mfj +(ImUy 



x 



d 3 k 1 



u e 



(13) 



which in the Boltzmann approximation reduces to: 



JR_ = Z^JULmTK (*t 
dM 2 ~ 7T 2 f 2 M 2 1 V T 



1 + 



2m 2 



1 



Ami 



i 

2\ 2 



M 2 J \ M 2 J 
ImU 

(M 2 -m 2 ) 2 + (ImUf 



(14) 



Following the notations of [3q we write the dilepton rate as: 



dR a 2 „„„ (M 
dM 2 = 2A^ MTK ^ ( T 



1 + 



2ml 
M 2 " 



(15) 



where we define 



F e// (M) ee 



m; 



48tt 

7 p 2 m 2 



/mil 



(M 2 -m 2 ) +(M) : 



(16) 
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To one loop[[!(J (i.e. consider only it — tt annihilation), we have 



"nn = -£Wl-±gV (17) 



Thus we recover the usual form of the form-factor, 



F eff (M)(l loop) = 1 



(M"-mlY + (m p r ptotal ) 2 V *PJ 
x9(M 2 -4ml) (18) 



where we've defined Imli = —m p T p tota i as usual. 

If we include other interactions besides n — n annihilation, in general 
I mil (or T p tota {) will increase, which increases the width of F e ff in Eq. (|T6D. 
However, the maximum of F e ff is given by 



48vrmJ 

F eff max - F eff (m p ) - - />2/mn (19) 

p p 



where we assumed that rh p and are slowly varying functions of M. Thus 
we see that the peak of F e ff decreases with increasing I mil. This contradicts 
the conclusions in |38| where increases in both the width and height of F e ff 



were reported. The discrepancy can partly be explained by double counting 
of dilepton rates in |39|, which is the basis for F e ff in [Q. To see this, we 



rewrite the dilepton rate fll3|) as 



dR 3 miT total 

(M)r^(M) /-d 3 k 



dM " 8 ^(M*-ml) 2 + (m p T pt otai) 2j » ^ 



(20) 
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where the decay widths are given byPTJP] 



r, 



i(M) 




/2/47T \3M 2 



1 - 



4m > 



2\ 2 



(21) 



L 2 m 2 



487T m p 
ImU 



m. 



1 - 



4m 2 



6 (M 2 - 4m 2 ) 



(22) 



We note that Q20D is essentially the same as the thermal p decay rate (2.9) 
of ||39|| , which was added to the 7r — 7r annihilation and other reactions in ||39|| . 
But Eq. ( P0| ) is simply a statement of the fact that the dilepton rate can be 
interpreted either as coming from annihilation of pions or as coming from 
the decay of thermal p's, the invariant mass distribution of which is given by 
the Breit-Wigner form (see also PD||[37||). We cannot, however, include both 
the thermal decay of p's and ir — ir annihilation. Doing so, as in [39 amounts 



to double-counting of the dilepton rates coming from the p channel. 



Note that T p ^ n7T in (2.18) of |37| has a factor of 16 instead of 48. 
In |4l[ this double-counting has been noted and corrected 
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